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Abstract 

The scattering approach to quantum transport through a disordered quasi- 
one-dimensional conductor in the insulating regime is discussed in terms of 
its transfer matrix T. A model of N one-dimensional wires which are coupled 
by random hopping matrix elements is compared with the transfer matrix 
model of Mello and Tomsovic. We derive and discuss the complete Fokker- 
Planck equation which describes the evolution of the probability distribution 
of TT^ with system length in the insulating regime. It is demonstrated that 
the eigenvalues of In TT^ have a multivariate Gaussian limiting probability 
distribution. The parameters of the distribution are expressed in terms of 
averages over the stationary distribution of the eigenvectors of TT^ . We 
compare the general form of the limiting distribution with results of random 
matrix theory and the Dorokhov-Mello-Pereyra-Kumar equation. 

PACS numbers: 72.10.Bg, 05.60,+w, 72.15.Rn, 73.50.Bk 
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I. INTRODUCTION 



The statistical properties of phase coherent quantum transport in mesoscopic systems 
have received increasing attention during the past few years [0. A variety of low tempera- 
ture transport quantities of conductors which are coupled to ideal leads can be expressed in 
terms of their scattering properties |2|-Q. Hence, their statistics may be studied in terms of 
the probability distribution of the scattering matrix. Several distinct approaches, including 
random matrix theory (RMT) ||-§|, Fokker-Planck (FP) equations [PHT3H, supersymmetry 
methods pl4|-|T6|| , and diagrammatic techniques p~7|JT8[] have been employed. This led to con- 



siderable progress in the understanding of quasi-one- dimensional wires fT9H23|| whose width 
is of the order of the mean free path which implies a structureless cross-section since no 
transverse diffusion takes place. The mean and the variance of the conductance of quasi- 
one-dimensional wires are now known for all length scales from the metallic to the localized 
regime |2"l-p6]. The generalization of these results beyond the quasi-one-dimensional regime 



is of considerable interest and has been the subject of some recent work |27|- |30|1 . Having this 



goal in mind, we focus on wires in the localized regime which are still quasi-one-dimensional 
in the sense that they are much longer than wide but which are not structureless in the 
transverse direction. In this regime the FP description simplifies considerably and progress 
is possible. 

This paper, which is the first of a series of two, deals mainly with the technical aspects of 
the problem and compares the general result which is obtained with previous results from 



random- matrix theory and the Dorokhov-Mello-Pereyra-Kumar (DMPK) [pl,10| equation. 
It has some overlap with the pioneering work of Dorokhov |32|] but goes beyond it by gen- 
eralizing the derivation of the transfer matrix limiting distribution for the one-dimensional 
wire by Kree and Schmid to the quasi-one-dimensional case. A preliminary account of 
the results which are presented here has been given in [34j|. In the second paper |35| we will 



investigate a model in which forward scattering is much stronger than backward scattering. 
We use the ratio of backward to forward scattering strength as a small expansion parameter 
and calculate the limiting distribution in the lowest two orders. 

The transfer matrix transforms the amplitudes of the propagating wave modes (open 
channels) at the Fermi energy in the left lead into the amplitudes of the right lead. A 
convenient parametrization for conductors with time-reversal invariance and with no spin- 



orbit scattering is the polar decomposition |32] , |36|| 




VTTx Vx 
Vx VTTx 




(i) 



where u, v are unitary N x N matrices and A is diagonal with real and positive diagonal 
elements Aj. The two-terminal conductance in units of e 2 / h is g = J2i % where % = 1/ (1+Aj) 
are the transmission eigenvalues of ttf and t = u(l/\/T + X)v is the transmission matrix. 

The transfer matrix of two samples which are joined together is the product of the 
transfer matrices of the individual samples. Building up a long wire by combining short 
samples, thus leads to a transfer matrix which is a product of a large number of random 
matrices. The eigenvalues of \n(TT^)/2L come in pairs (a m (L), —a m (L)) where 1 + 2A m = 
cosh(2a m L) = coshr m and L is the system length. The corresponding eigenvectors are 



y u rm a rr, 



/\/2 and (u^, —u* l J) T /\/2 where u m is the m th column vector of u. 
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From Oseledec's theorem |37| for random matrix products it is known that the a m (L) 
are self- averaging and distinct if L goes to infinity. The limiting values = lim^oo a m (L) 
are the Lyapunov exponents. They characterize the rate of exponential growth of the A m 
with system length. Furthermore there are central limit theorems [0,0 which show that u 
and v have stationary distributions and that the quantities (r m — 2a^L) / \[L have Gaussian 
limiting distributions if L goes to infinity. Oseledec's theorem implies that the T m can be 
ordered as 1 < T 1 < T 2 ■ ■ ■ < T N if 2afL > 1. In this regime g = £ m 2/(2 + cosh T m ) pa 
4exp(— Ti) = 4exp(— 2a\V) and the sample is strongly insulating. The decay length £ = 
l/2a^° of the typical conductance is usually identified with the localization length. Johnston 
and Kunz [^(J applied the central limit theorems to the Anderson model. We derive the 
Gaussian limiting distribution within a FP approach, which will establish a link between 
the parameters of the limiting distribution and the stationary distribution of u. 

The determination of the Lyapunov exponents of random matrix products is a problem 
which arises often in the context of disordered systems. At present there is no powerful 
method to calculate them analytically. Only special cases like sparse random matrices 



have been solved ||41|| . The weak disorder expansions which have been developped cannot be 
succesfully applied to quasi-one-dimensional conductors because of the problem of degenerate 
eigenvalues P2}415[]. The full limiting distribution has been mainly studied in numerical 



simulations [46,47|. Apart from the numerical data there are only two analytical approaches 



which make quantitative predictions, RMT p|,|6|,^8| and the DMPK equation |22|.|2~ 



The RMT ensemble maximizes the information entropy of the transfer matrix probability 
distribution subject to the constraint of a given density (p(T)) = (X) m 5(r ~ T m )). As 
a consequence, T and the unitary matrices are stochastically independent. The unitary 
matrices are isotropically distributed and the probability distribution of T has the form 

p(T)=X- 1 ex P {-PH(T)}, (2) 

where 

H(T) = - ^ | cosh T m - cosh r n | - ^ £ ln(sinh T m ) + £ V(T m ) (3) 

m<n r m m 

and Af is a normalization factor. The parameter (3 is determined by the symmetry of the 
transfer matrix ensemble. The orthogonal {(3 = 1), unitary ((3 = 2), and symplectic (f3 = 
4) ensembles correpond to conductors with time reversal symmetry, broken time reversal 
symmetry, and strong spin-orbit scattering respectively. The potential V(A) has to be 
determined from (p(T)). RMT is known to be a good but not exact description of quasi- 
one-dimensional conductors without transverse structure 0. For such conductors and large 
N, (p(T)) ta NI/2L if < V < 2L/1 and (p(T)} ta if T > 2L/1 where / is the mean free 
path. The resulting potential is quadratic V(T) ~ NIT 2 /AL [[49]. In the insulating regime 
where 1 <C Y\ C T 2 --- C Tjv the Hamiltonian (0) simplifies, since ln(sinhr m ) ps Y m and 
In | coshr m — coshr„ | pa T n if m < n. This leads to the Gaussian probability distribution 

where (T m ) = (m - 1 + 1//3)2L/(IN) and a 2 = 2L/(/3lN). Note that any potential must 
have the form V(r) ~ aT + bT 2 if L ^> £ and T ^> 1 in order to be consistent with 
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the Gaussian limiting distribution. This implies always equidistant mean values T m and 
equal variances for the fluctuations around them. Numerical simulations of conductors with 
transverse structure show that this is in general not true |5(],|46|.|4^J . Therefore RMT cannot 
describe such conductors. 
The DMPK equation 

M£E> « (|L +fl *SE)), (5) 

oL ll ^xoYm \dY m dT m J 

where fi(T) = — J2 m <n h 1 |(cosh Y m — cosh T n )/2| — 1/(3 J2 m In | sinh T m | and 7 = (3N + 2 — (3, 
constitutes an exact description of quasi-one-dimensional wires without transverse structure. 
Its solution 

P (T)(x n |coshr m -coshr n |^ 2 |r^-r^| 

m<n 

x n[exp(-70/8L)r m (sinh r m ) 1/2 ] (6) 



in the insulating regime p2| , p3| can be as well approximated by a Gaussian distribution of 
the form (|) if 1< V x < Y 2 ■ ■ • < V N} where (T m ) = [1 + (3{m - 1)}2L/[1((3N + 2-/3)] and 
a 2 = AL/[l((3N + 2 - f3)}. Note that the mean values of RMT and of the DMPK equation 
coincide for large N whereas the variances differ by a factor of two. 

The content of the paper is organized as follows. In Sec. [TI| a Hamiltonian model of N 
one-dimensional wires which are coupled by random hopping matrix elements is compared 
to the transfer matrix model of Mello and Tomsovic [[27],[2l|. In Sec. I J I] we derive the 
FP equation which describes the evolution of the probability distribution of T and u with 
system length in the localized regime. In Sec. [TV] we generalize the derivation for the transfer 
matrix limiting distribution of a one-dimensional wire by Kree and Schmid |33[] to the quasi- 
one-dimensional wire. A first application of this approach is presented in Sec. |V], where 
we investigate the Equivalent Channel Model (ECM) of Mello and Tomsovic. The joint 
probability distribution p(L; T) of this model is known to be identical to the distribution of 
the DMPK equation for (3 = 1. We recover the Gaussian distribution ([|) and show that the 
stationary distribution of u is isotropic. 

There are four appendices. The derivation of the FP equation in Sec. [I I J is based on 
a simplified version of the general Langevin equations for Y and u which are obtained in 
Appendix [AJ The coefficients of the FP operator are derived in Appendix |B[ In Appendix 
we show that a measure for the unitary group which has been introduced in the text is the 
invariant measure. An alternative derivation of the FP equation is presented in Appendix 
0. The summation convention is used throughout the whole paper. 



II. HAMILTONIAN VERSUS TRANSFER MATRIX MODELS 

The FP approach to disordered conductors has been pioneered by Dorokhov | 32| . He 



started from a microscopic model of disordered coupled chains which led to a quite compli- 
cated FP equation. Similar models were also studied by other techniques PH-p^ . Recently 
Mello and Tomsovic proposed a class of models which was formulated on the level of the 
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transfer matrix p7,28|. On the one hand these models lead to simpler FP equations but 



on the other hand the underlying Hamiltonian is not known. In this section we propose 
a microscopic model which is simpler than the one which has been used by Dorokhov and 
compare it to the model class of Mello and Tomsovic. 

Consider the scattering of electrons at a quasi-lD disordered conductor with a (d — 1)- 
dimensional cross section which is connected to perfectly ordered leads. The conductor 
consists of N ID-wires which are only coupled by random hopping matrix elements. It is 
described by the Hamiltonian 

h 2 

H nn > = -<Wo 9 2 x + V nn '{x), (7) 

where V nn i(x) is real and symmetric in its indices and n = 1, ■ ■ ■ , N. The potential V nn i(x) 
is zero in the leads and stochastic in the disordered system of length L. It describes on-wire 
disorder for n = n' and random hopping between the wires for n ^ n'. The independent 
matrix elements of V(x) are chosen to be uncorrelated and Gaussian distributed with zero 
average 

(V rm >{x)) = 

{V nn >(x)V mm >(x')) = U nn '5(x - x')(6 nm 6 n i m > $nm'$n'm) (8) 

where U nn > = U n > n . The special case U nn i = U/N can be interpreted as a continuous one- 
dimensional iV-orbital model [|56|1 , which is connected to ideal leads with no exponential 
decaying modes. 

The solution *ff(xn; E) of the scattering problem with the incoming waves 



¥ n {xn- E) = ^Jjj- (a l n exp{ikx} + b r n exp{-ikx}) (9) 

is an eigenfunction of the Schrodinger equation with energy E = h 2 k 2 /2mo. Its form in the 
left and the right lead respectively is 



¥ /r (xn; E) = (af exp{^x} + hf exp{-ikx}) (10) 

where the amplitudes a|/ r and b 1 ^ have been normalized in such a way that the probability 
current in x-direction is j x = Y^n \ a n\ 2 — \b n \ 2 - The S-matrix transforms the amplitudes of 
the incident waves into the amplitudes of the scattered waves 

(11) 

Current conservation and time reversal invariance imply that S is unitary and symmetric 
respectively. The transfer matrix by contrast transforms the amplitudes in the left lead into 
the amplitudes in the right lead 








t-r't'^r r't'- 1 
-t ' V t'- 1 



(12) 
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Here, current conservation and time reversal invariance leads to the form 




(13) 

where oca" 1 — f3/3' = 1 PJT0[1. Apart from the polar decomposition (fl|) there is another useful 
parametrization of the transfer matrix which has been introduced by Mello and Tomsovic 
and has the form 



exptf \ / y/1 + 7777* 77 
exp i9* / I 77* y/1 + 77*77 



(14) 



where t3 and 77 are complex NxN matrices and = — and 77 T = 77. The wave amplitudes 
ai/ r and bi/ r fix the values of \l/(x, n; E) and d x ^(x, n; E) at the edges of the sample. This 
implies that the transfer matrix of two samples which are matched continuously together is 

T(L + L',0)=T(L',L)T(L,0). (15) 

Hence, the transfer matrix of a sample of length L can be obtained by dividing it into short 
segments of length 5L and multiplying the transfer matrices of the segments. The evolution 
of the transfer matrix with the system length is a multiplicative stochastic process. It can be 
described by a Langevin equation since the model is continuous in the scattering direction. 
The Langevin equation has the form 



dT(x, 0) _ f ~f u {x) 7 12 (x) 
dx ~ [ l 2 \x) 7 22 (z) 



T(x,0) (16) 



with the noise 7 u (x). The symmetries 

f 

7 21 =7 12 *, (17) 



22 11 * 



and 



7 Ut = -7 n 

,12 T 12 



7 12T = 7 12 (18) 

which will be derived below ensure time reversal invariance and flux conservation. Iterative 
integration of the Langevin equation yields 

T(x + 5L,x ) = 1+ ( 21 £22 (19) 



where 



s ij = dx~r lj (x)+ dx dx'Y k (x)-r kj (x') + ---. (20) 

J XQ Jxq J Xq 
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For uncorrelated noise the first term of this expansion is of order [SL) 1 / 2 and the second term 
is of order SL. For the derivation of the symmetries flTTj ) and flTBp, however, it is convenient 
to start with a finite correlation length of the noise. Taking Sx < SL and SL to be smaller 
than this correlation length one may expand ^(xo + Sx) = ^(xo) + d x -f^ (xo)Sx + 0((Sx) 2 ) 
which leads to 

s ij = ^{x )SL + O{{SL) 2 ). (21) 

Equation (|T3|) enforces the symmetries 7 22 = and 7 21 = 7 12 *. Comparing the expan- 
sion ([H]) with the parametrization (|T4]) of the transfer matrix one finds 

■&(x + SL, x ) = 7 n (x )5L + 0(5L 2 ) 

r](x + SL, x ) = j 12 (x )SL + 0{SL 2 ) (22) 

which implies 'y 11 ^ = — 7 11 and *y 12T = -y 12 . These symmetries remain valid in the limit of 
zero correlation length. 

In the sequel we derive 7 n (xo) an d 7 12 (^o) f° r the Hamiltonian model (0). The stationary 
solution of the scattering problem obeys the Lippmann-Schwinger equation 

V{xn;E) = m m (xn-E) + / dx x £ G + (m,x ini ;£)K in ;( 3 ; 1 )1'(a; 1 ,n' 1 ;£), (23) 



m,n' 



where 

G+(xn,x'n';E) = e^{ik{E)\x - x'\] (24) 

is the free retarded Green's function. Iteration of the Lippmann-Schwinger yields the Born 
series 

$(xn; E) = ¥ n (xn;E) 

+ / d Xl G+{xn,x 1 n 1 ;E)V nin ,{x 1 )W n {x 1 ,n' 1 ;E) + --- (25) 

which can be translated into series for the transmission and reflection matrices by Eq. ([|) 
and Eq. (pTTJ) . The n th orders of these series are at least of the order (SL) n since they contain 
n integrations from Xq to Xq + SL. Thus, only the first orders can contribute to the terms 
of order SL of the expansions t = 1 + fiSL + ■ • ■ , r = r Y SL + • • • , t' — 1 + t' l SL + • ■ ■ , 
r' = r' 1 SL + • • • . Inserting these contributions into the relations 7 11 = t 1 and 7 12 = r' 1 
which follow from Eq. ( ]T2| ) and Eq. fl2~T|) yields 

11 / \ -im 

7nn'(^o) = -TsrKn^o) 

12 / \ -im exp(-i2kx ) , , , , 

Tnn'W = T2T V nn ix ). (26) 

The phase exp(— i2kxo) which appears in 7 12 (xo) is a consequence of the transformation 
rule 
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T(L + x ,x ) 



e -ikx 01 Q \ / Q ikx 01 







(27) 



which accounts for a shift of the disordered region by xq. 

Now we are in the position to compare the Hamiltonian model with the transfer 

They divided a sample of length L into 



matrix model of Mello and Tomsovic [ pY| , pq 
n = L/5L uncorrected scattering units with identical statistical properties. Specifying the 



first two moments of $ m „ and rj mn for one scatterer 



5L = (f]mr 


)SL 


= 


(ftmn&m'n 


ISL 


11,11 

— K 

mn,m 


{$mn' l 9 m 'n 


)<5L 


11,22 
^mn,m 


(jlmrJlm 1 n 


ISL 


12,12 

— K 

mn,m 


(VmnVm'n 


)SL 


12,21 

— K 

mn,m 


($mnVm'n 


)5L 


11,12 

— K 

mn,m 


($mnVm'n 


)SL 


11,21 

— K 

mn.m 



(28) 



and taking the continuum limit of a high number of weak scattering units such that 



q ^J^ mn,m'n' u mn,m'n' 



5L 



(29) 



and that 1/SL times higher moments gives zero in the same limit determined completely 
the stochastic evolution of the transfer matrix. As a consequence one finds 



11 1 

mill 



12 1 

rnnl 



( 11,11 4- 12 ' 21 "l II 
v^mk.kn ' ^mk,kn) I 



a 



.11,12 

mk,kn 



13 F i>j ' 1 



ijA'j' 
mn,m'n' 



(30) 



where [• ■ ■] = ]imsL->o(- • -)sl/^L. The same limit for higher moments of is zero. Mello 
and Tomsovic have chosen the following simple model for one scattering unit. Assume that 
the independent matrix elements of i? and r/ are uncorrelated and that their phases are 
randomly distributed. Averaging over the arbitrary distribution of their modulus then leads 
to 



\e ij 1 



If 11 f 11 1 
L mn~m'n'\ 



r 11 22 ■ 

I mn m'n' . 



\f 12 f 21 1 
L mn°m'n' 1 



1 1 

"mn 1 "nm' f 

Imn 

r c J_ 

"mm' "nn' r 

Imn 



Omm 1 Onn' ~t~ O mn i 5 nm i 1 



l + 6„ 



It 
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If 11 f iz = n 

L c mn c m'n'J 

[c« F 21 I = 
\r mn^ ra'n' \ w 

[ £ mn £ m'n'] = (31) 



where and l b mn are the mean free paths for forward and backward scattering from channel 
m into channel n and l/l^ b = Yin lAmn are the total inverse mean free paths. The inverse 
mean free paths 1/P mn and l/l b mn are defined by the probabilities per length [|i mn — <5 mri | 2 ] and 



r 



mv | 



2 ] for a forward and a backward scattering process respectively. By Eq. (^) and Eq. 
(0) [|*mn — <^mn| 2 ] = [| £ mnl 2 ] an d [l r mn.| 2 ] = [| £ mJ 2 ] wmcn leads to the above identification 
of the model parameters with the mean free paths. 

Now we calculate [e^J and \^i in e% m'n\ f° r the Hamiltonian model (§). Inserting 7 11 (x ) 
and 7 12 (xo) from Eq. (P5|) into Eq. ( |2"0| ) and averaging over the Gaussian white noise yields 

[e ij 1 = 

T ij i'j' 1 ^i^jj ^ J ) ^mm'^nn' ^mn'^nm' ( r >0\ 

[ £ mn £ m'n'\ = / i , r l^J 

where l/7 mn = l/l f mn = l/l b mn = (m /{h 2 k)) 2 U mn (l + 5 mn ) and 1/7 = Y m l / l mn- Note that 
V mn and are not independent as in Eq. ([31]). The coefficients c(ij.i'f) which are not 
related through the symmetries of ^{xq) are c(ll, 11) = — c(ll,22) = — c(12,21) = — 1, 
c(ll,21) = — c(ll, 12)* = exp(i2kx ), c(12, 12) = — exp(— i4kx ). Hence, the moments 
[ £ mn £ m'n'] which vanish in model (|3T|) oscillate with x in Eq. (|32|). This will cause the 



coefficients of the FP equation (^) for the probability distribution of TT^ to oscillate with 
the system length. In the limit of weak disorder {kl 1) these oscillations are very fast on 
the scale of the mean free path I which is the characteristic length over which the probability 
distribution changes. Then, it is justified to average over the oscillations, which amounts 
to replace the oscillating moments in Eq. ([32]) by zero. The resulting model is very similar 
to the model fl3"T|) if P mn = l b mn but not equivalent. It would be equivalent if the phases of 
Vmn were not random but had been chosen to take the values exp(«7r/2) and exp(— m/2) 
with equal probability. For stronger disorder it is no longer justified to average over the 
oscillations. 

As an alternative to the continuum limit of Mello and Tomsovic one may specify directly 
the statistics of Ymn( x )- Choosing Gaussian white noise such that 

(lmn( X )lm'n'( X ')) = ^mn!m'n'( X ) S ( X - x ) (33) 

&0o)] = <r%%(x )/2 

[ £ mn( X o) £ m'n'( X o)} = °%n?m> n> : o) • (34) 



leads to 



The Hamiltonian model (g) and the model ( pT|) are special cases of this class of models. 
Note however that [e^ n (x )] differs in general from the result ( |3"0p of the continuum limit. 
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III. LANGEVIN AND FOKKER-PLANCK EQUATIONS 



The evolution of the transfer matrix with the system length is a stochastic process 



which can be described by Langevin and FP equations. Dorokhov recognized that the 
stochastic equations for the matrix 



M = TT ] 




( coshT sinhT 
I sinh r cosh T 




(35) 



are closed which allows to eliminate the degrees of freedom of v. We follow Dorokhov and 
start from his Langevin equations ( |A10p for T and u which are derived in Appendix |^ 



for the sake of completeness. Due to the self-averaging of the a m one expects that in the 
insulating regime where L ^> £ = l/(2a^°) the r m can be ordered: 1 <C 1^ <C ■ ■ ■ <C T N . 
This ordering justifies the neglect of exponentially small contributions to the terms 

cothTj = 1 +0(exp(-2r i )) 

sinhr n = f + 0(exp(r n - Tj)) n<j 
cosh r n - cosh Tj { 1 + olexplTj - T n )) n>j 1 ' 



of the general Langevin equations ( A10 ) which leads to the considerable simplification 

E mm -\- E 



dT 



dL 

du 
dL 



7 n w + uP. (37) 
The matrix elements of P are P mn = 6(n — m)E mn — 9{m — n)E mn where 




e( n ) = { \ n = (38) 



and E = u^ 12 u*. The symmetries 7^ = —7 and P^ = —P imply that d(uu^)/dL\ uu ^ =1 = 
0. Thus, the simplified Langevin equations (|37D still conserve the unitarity of u. The 
stochastic process described by them leads to a limiting distribution for (r m — 2a^L) / y/T 
which is independent of the initial conditions. Hence, this limiting distribution must be 



identical to the one which is produced by the original Langevin equations (|A10|) as long 
as 1 T 1 ■ • • <^ T N . Therefore, it is possible to use the simplified Langevin equations 
together with convenient initial conditions to determine the form of the limiting distribution 
in this parameter range. 

Due to the neglect of exponential small terms in fl36|) , r m can become negative and it 
is natural to extend the range of T m to —00 which is justified because the probability to 
find a negative value of T m will turn out to be exponentially small. For similar reasons we 
relax the strict ordering of the T m . A parametrization of it by a set of N 2 independent 



10 



parameters seems to be rather complicated. Instead, we extend the range of the matrix 
elements u mn = x mn + iy mn to arbitrary complex numbers thereby obtaining a stochastic 
process on a higher dimensional cartesian space. The standard derivation technique [[57]] for 
the FP equation of such a process yields 

dp(L;T,u,u*) = {dTmdr J mn + dvm z m + d)p{L . r u (3g) 

where l\ m dT m Y\ m i n i dx m i n idy m i n i is the measure of the cartesian space. The operators A mn , 
B m , and C are 

A m n ^t^"^" 1 ^"^™] 

B m = -[AT m ] + d Um , n , [AT m Au mW ] + d u * , , [Ar m A 

lib It J m'n' 



U 



m'n' 



C = -d Umn [Au mn ] - d u * mn [Au* mn ] + d Umn d u * [Au mn Au* m , n ,] 

mil m'n' 

+\d Umn d Um , n , lAu mn Au m , n ,} + )^d u * mn d u * mini [Au* mn Au* m , n ,}, (40) 

where Ar m = T(L + 5L) — T(L) and Au mn = u mn (L + 5L) — u mn (L). The brackets [■ • ■] 
define the coefficients of the FP operator and stand for \imsL-*o(- ■ -)sl/^L where (• • is 
the average over the disorder in the region between L and L + 5L. The explicit form of the 
coefficients is derived in appendix |B[ 

The multiplication of a probability distribution on the cartesian space by the 5-function 



) n *(i>( )Wem 



u m'nU n i n ) 



m=l \ n / rn'<n' \ n 



(41) 



restricts it to the unitary group. Since the Langevin equations (|37|) conserve unitarity, one ex- 
pects that the FP operator commutes with the ^-function. In fact, the operators A, B and C 
commute with every function of the type f{uv)) for arbitrary complex matrices u. A lengthy 
but straightforward calculation with the coefficients ( |B2| ) which exploits the symmetries 
(|I~8l) proves that this is true. Thus, the restriction to the unitary group may be incorporated 
into the new measure rim dT m dfi(u) where dfi(u) = V _1 (A^)5(1, uu*) Y\ m i n i dx m > n >dy m ' n > and 
V(N) is the volume of the unitary group. It is shown in appendix that dfj,(u) is the invari- 
ant measure of U(N). We note that an integral of the type / dfi(u)Xg(u, u*) (X = A, B, 
or C) with respect to the invariant measure is evaluated in two steps. First, the 5-function 
is commuted with X. Second, the integrations are carried out yielding only contributions 
for terms which do not have derivatives with respect to u mn or u* mn in front of them. 

It is worth emphasizing that the operators A, B and C do not depend on T. This is 
the great simplification which has been achieved by the neglect of the exponential small 
terms in (0). However, they can still depend on the system length L as it is the case for 
the Hamiltonian model (0). There, the factor exp(— ikxo) of 7 12 (^o) in Eq. ( p6|) leads to 
terms which oscillate with the system length. For weak disorder (kl ^> 1) the oscillations 
are fast on the scale of the mean free path / which justifies replacing the oscillating terms by 
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their averages. For stronger disorder the oscillations can be absorbed into the new variable 
Umn = u mn exp(ikL). The transformation of the FP equation to this variable leads to the 
additional term —ik{u mn du mn — u* mn du* mn ) in the FP operator. The resulting FP equation 
is very similar to the FP equation of Kree and Schmid [ 33f] . It is the generalization from 
N — 1 to arbitrary channel numbers. 



IV. THE LIMITING DISTRIBUTION OF THE TRANSFER MATRIX 

Kree and Schmid discussed thoroughly the asymptotic probability distribution for the 
Landauer conductance g = \t\ 2 of a long one-dimensional wire The transmission and 
reflection amplitude for incident waves from the right are t = u(2/(l + coshr)) 1 / 2 ^ and 
r = u 2 ((coshr — l)/(coshr + l)) 1 / 2 where u and v are simply phases. They showed that 
(T — 2a°°L)/^/L has a Gaussian limiting distribution if L — > oo which implies a log-normal 
distribution for the conductance. The parameters of the Gaussian distribution were ex- 
pressed in terms of averages over the stationary distribution of u. Similarly we expect 
that the corresponding quantity (T — 2a oc L)/\/ r L for the quasi-one dimensional wire has a 
multivariate Gaussian limiting distribution whose parameters can be expressed in terms of 
averages over the stationary distribution of u. 

It is useful to look at the first moments (r m ) l in some detail before deriving the general 
form of the limiting distribution. Integrating Eq. fl3~9|) with respect to T leads to the closed 
FP equation 

dq (L;u,u>) =dq(L . uu , ) (42) 

for u. We expect that the stationary distribution of u is the unique stationary solution 
Q s tat{u,u*) of this equation. Hence, the spectrum of C should consist of one eigenvalue vq 
which is zero and others z/, with negative real parts. The smallest absolute value of these 
real parts is henceforth called v. In the sequel it will be assumed that the eigenfunctions 
qi(u, u*) of C form a complete system where qo(u, u*) = q s t a t(u, u*). Then, any probability 
distribution may be expanded into X)i Ciqi(u, u*) where conservation of probability implies 
/ dfj,(u)qi(u, u*) = for i ^ and Co = 1. Multiplying Eq. with T m and integration by 
parts yields 

d ^Q T f L = ~ J dfi{u)B m q(L; u, u*) 

= J dfi(u)[Ar m ]q(L;u,u*). (43) 

Expanding the initial distribution into eigenfunctions q(0;u, u*) = q sta t(u,u*) + X^o c « 
qi(u, u*) gives q(L; u, u*) = q s tat(u, u*) + J2i^o <k exp(uiL)qi(u, u*). This leads to the large 
length asymptotic behaviour 

(T m ) L - (r m } w L J dfi(u) [AT m }q stat (u, u*) + [const. + 0(exp{-uL})} (44) 

where the terms in brackets result from the crossover of the initial into the stationary distri- 
bution of u. The self-averaging = limz / ^ 00 a m (L) = limL^ 00 (a m (L))i of the Lyapunov 
exponents implies = \im L ^ QO (T m ) L /2L which leads to 
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-J d/i(u)[Ar m }q stat (u,u*). (45) 



This relation has been first derived by Dorokhov | 3~2|| . Using Eq. ( |B2|) to calculate [Ar r 
for the model (|3l|) yields 

9(m — k 3 ) f / 

a ™ = Wb Trr* V / "M") [ U kimU k 2k 3 U klk3 U k2m + U kim U k2k U klk3 U k2m (46) 

^fcifcal 1 + d kik 2 ) J 



+^u kim u k2k3 u* kim ul 2k \ q sta t(u, U*). 



The same formula has been obtained by Chalker and Bernhardt |29j for the special case that 
there is only back scattering into the same channel. They discussed also the consequences 
of this relation in the context of the Anderson transition. 

In the sequel we will go beyond the first moment (T m ) l and derive the general form of 
the limiting distribution of (T — 2a°°L) / \J~L. For the sake of simplicity we choose the initial 
distribution 

N 

p(0; T, u, u*) = n S{T m )q sm {u, u*) (47) 

m=l 

which implies that (T m ) L = 2a"^L (see Eq. (JUD). The formalism which will be devlopped 
below could be used to show that different initial condition would not change the form of 
the limiting distribution but only the way it is approached. It is convenient to introduce 
Dirac notation 

(01(7=0 (7|0)=0 
r=\0)(0\ Q=l-V (48) 

where q sta t(u,u*) = (u,u*\0) and (0\u,u*) = (0 1 0) = 1 so that the average J dfi(u) 
Xq stat (u,u*) may be simply expressed as (0|X|0). 

The central quantity which will be used below to derive the limiting distribution is the 
generating function 

P(L; r, u, O = I J] dT m , exp{z ]T r m {T m - (T m ) L )}p(L; T, u, u*) (49) 
m' m 

for the central moments 

C L {miri, • ■ • , m k r k ) = ({T mi - (T mi ) L ) ri ■ ■ ■ (T mk - (T mk ) L ) rk ) L 

= J ^{u){-id Tmi r--- (-zd Tmk Y*P{L- r, u, «*) | T=0 . (50) 

The Fourier transform of the generating function gives back the probability distribution 

i r N i N \ 

p(L;T,u,u*) = / ]J dT m exp<-iJ2( r n- (r n )L)T n > P(L;t,u,u*). (51) 

\ Z7V ) J m=l { n=l ) 

This implies the evolution equation 
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dP(L;T,u,u ._ r> „ 



<9L 

where 



~~ r m T nA mn — iT m B m + C)P(L ;t,u,u*) (52) 



d(T m ) L 



dL 

B m - (0\B m \0) (53) 



and P(0; r, u*) = q s tat{u, u*). The formal solution of Eq. (52) is 

P(L] t, u, u) = exp{(-T m r n A mn - iT m B° m + C)L}q stat (u, u*). (54) 

We follow Kree and Schmid and reexpress it in terms of the operator generalization of the 
Cauchy formula exp(5L) = l/(2m) § d( exp(i(L) / (( + ib), 

P(L] t, u, u*) = j> d( exp{i(L}R(r, Qq stat (u, u*) (55) 

where 

R(t, C) = (C + i{-T m T n A mn - iT m B° m + C))' 1 (56) 

is the resolvent operator and the integration contour encircles all the eigenvalues of 
—i(—T m T n A mn — ir m B° m + C) counterclockwise. Equation (^) for the moments of (r m — 
(r m )_L) only requires the resolvent operator for infinitesimal values of r. The spectrum of 
— £(— T m r n A mn — iTmB^ + C) then lies in the neighbourhood of the spectrum of —iC which 
consists of one zero eigenvalue and other eigenvalues in the upper half plane. Hence one can 
choose the integration to run just below the real line from minus to plus infinity and close 
the contour in the upper half plane. 

We do not reconstruct the limiting distribution from the moments C^miri, • • • ,mkTk) 
but proceed indirectly. First consider the linear combination Q = c m Q m = c m (T m — (r m )i) 
with arbitrary coefficients c m and study its moments which are given by 

{Vl r } = J d^(u)(-id T ) r P(L;rc,u,u*)\ T=0 

= ^ f d( e xp{t(L}(-td T y(0\R(rcX)\0)\r=o- (57) 

R(tc, C) = (C + t(-r 2 c m c n A mn - iTc m B° m + C)' 1 (58) 



Since 



is similar to the resolvent operator of a one-dimensional wire one can employ the technique 
of Kree and Schmid to calculate the moments. Expanding R(tc, () into powers of r yields 

oo . 

R(tc, C) = i? (C) E [(*T 2 c m c n A mn - Tc m B°jR {()] ' (59) 

fc=0 
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where -Ro(C) = (C + iC) 1 rnay be decomposed into a part which is singular at ( = and a 
non-singular part 

Ro — Ros + Ron 

= - + ~^-r. (60) 

( ( + iC 

Only the terms of order r r of the expansion (59) contribute to the r th moment of Q. The 
first moment 

(O)i = /dCexp{iC^}(0|A)(C)ic^mA)(O|0) 

Z7T2 J 

= ^(0K|0)/dC^^i (61) 



is zero because of Eq. (|53"D as it should be. Collecting the terms which contribute to the 
second moment yields 

(^ 2 ) = ^^dCexp{iCL}(-2c m c n )(0|iA)>lmnA) + A)^A)^°A)|0) 

= ^/rfCexp{zCL}^^(0|a mn + 2&_^i£|0) (62) 

The residues of the pole of second order at ( = are 2c m c n (0\A mn \0)L and —2c m c n 
(tyB^C^ 1 B^L + £?^C ,_2 -B°|0) for the first and the second term respectively. The poles in 
the upper half plane of the second term only give rise to exponential decaying contributions. 
Hence 

(n 2 ) = 2uL + const. + O (exp{-z/L}) (63) 

where 

u = c m c n (0\A mn -B° m C- 1 B° n \0) 

C-mC-n-Amn- (64) 

It can be shown along the same lines as in work of Kree and Schmid that the higher moments 
have the form 

n\ 

(Q 2n +i) = 0{L n ). (65) 

For details of the proof we refer to [Q. The form of the moments implies that VL = Q/yL 
has the Gaussian limiting distribution 

Soo (n) = -^exp{-{l 2 /(4Lu)}. (66) 
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Now we discard the finite length corrections to the limiting distribution and keep only 
the universal parts (n 2n )^ = (ujL) n (2n)\/n\ and (Vt 2n+1 )^ = of the moments (|65|). The 
corresponding universal part of the generating function S(L; r) = / dfi(u)P(L; r, u, u*) 
is denoted by S u (L;t). Since (Vt r ) u L = (-id T ) r S u (L;rc)\ T=0 we find S u (L;rc) = 
exp{— T 2 c m c n A m nL} which implies 

S U (L; t) = exp{-r m r n A mn L}. (67) 

The Fourier transform (|5T| ) of S U (L; r) yields the universal part of the probability distribu- 
tion 

s u (L;T) = - — 1 exp{-(r m - (T m ) L )A^ n (T n - (T n ) L )/AL} (68) 

(47rL)^/ydet({An4) 



where (r m ) i = 2a°°L. Hence, the limiting distribution of fl = (T — 2<x co L) / \fL is 

Soo (fl) = / exp {-n m A m ln n /4} . (69) 



Note, that the form of the correlator 

((a m - OK - 0>l = -^{-id Tm ){-id Tn )S u {L- r)| r=0 

~^\Aran Anm) (70) 

implies that the a m = T m /2L are self- averaging and that the fluctuations around their 
limiting values are in general correlated. Such correlations are not predicted by RMT or the 
DMPK equation but are consistent with numerical simulations p6|,|47| . 

The variance of the Gaussian distribution of In g « — Fx + const, in the insulating regime 
which follows from var(lng) « var(I\) and Eq. ( j63|) is 

var(ln^) w 2A n L + 0(l). (71) 



V. THE EQUIVALENT CHANNEL MODEL AS A SPECIAL CASE 



Mello and Tomosovic have shown [ [ZT| , [ZB| that the joint probability distribution s(L;T) 

(72) 



of the ECM which is the model (|31|) with backscattering mean free paths of the form 

1 1 + 5, 



obeys the DMPK equation for (3 = 1. The form (^) of the solution for the DMPK equation 
in the insulating regime is a special case of the multivariate Gaussian distribution (|68|). 
Hence, we expect to recover this solution from our approach if we apply it to the ECM. 

Evaluating Eq. (|B2|) for the coefficients [Ar m Ar n ] and [Ar m ] with the backscattering 
mean free paths (|72|) yields 
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2 

[Ar m Ar„] = JfoTjy _|_ ^\ {y , k 1 m U k 2 rn U kin u k2n + u kim u k2m u k 1 n U k 2 n 



[ AT m] = (Kim U *k2k 3 u kxk 3 U k2 m + Wfc 1 m«fc 2 fc 3 w fc 1 fc 3 M tem ( 73 ) 

~^~ U kim U k2k 3 U k2k3 U kim "I" u k\m u k2k2, u k2k 3 U k\r 



Since the coefficients have the constant values [Ar m Ar n ] = 45 mn /(l b (N + 1)) and [Ar m ] = 
2m/ (/ b (iV+l)) if u is unitary, one can integrate the evolution equation (|52|) for the generating 
function P(L; t, u, it*) with respect to w. The solution of the resulting equation 



OS{L-t) 2S mn 

-^ = -nNTT) TmTnS{L]T) (74) 

is S(L; t) = exp{—2J2m T mL/(l b (N + 1))} which leads to the Gaussian distribution 

s(i;r) = n^-pp^ 2 ) (75) 

with (r m )i = 2mL/(l\N + 1)) and ct 2 = 4L/(l b (N + 1)) as expected. 

It is a specific property of ECMs that the evolution of the joint distribution s(L; T) 
decouples from u. Therefore, the limiting distribution does not depend on the station- 
ary distribution of u. Still, it is of interest to know the stationary distribution. Solv- 
ing the equation Cq sta t{ u ,u*) = which determines it is in general a difficult problem. 
We demonstrate below that due to the simple form of the backscattering mean free paths 
the stationary probability measure q s t a t{u, u*) flm,™ dRe(u mn )dIm(u mn ) can be found to be 
V _1 5(l, uu') Ilm.n dRe{u mn )dIm{u mn ) which is the invariant measure of the unitary group. 
Using the form ( piO[ ) of C in which the derivatives act directly on the distribution and the 
property that C commutes with the 5-function yields 

(75(1, uv?) = <f(l,tmt) Qmn[Ar m Ar n ] - m[AT m ]) 

= 5{l,uu ] )g(uu ] ) (76) 

where g(uu^) is zero if u is unitary. Applying higher powers of C to the ^-function gives 
C k 5(l, wu)) = 5(l,uu^)g(uvJ) k . Hence, the initial distribution q(0;u,u*) = V _1 5(l, uw') 
evolves into 

q(L; u, u*) = V _1 5(l, W) exp{g(uu^)L} 

= V -1 £(l, tm + ) (77) 

which shows that it is the stationary distribution. 



VI. CONCLUSION 

The general form (^) of the limiting distribution and the link between its parameters 
and the stationary distribution of u are the main results of this paper. This form implies 
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that the RMT probability distribution (||) is not sufficient to describe quasi-one-dimensional 
conductors with transverse structure. The generalization of RMT to such conductors remains 
a challenging problem ||58|| . Beenakker [22] has shown that a correct description of quasi-one- 
dimensional wires without transverse structure requires a modification of the interaction in 
the Hamiltonian (^). It is not clear whether a modification of the interaction is sufficient to 
describe conductors with transverse structure or if three and more eigenvalue interactions are 
needed. Since any generalization must be consistent with the form (^) it is of considerable 
interest to have explicit results for (T m ) and A mn which go beyond RMT and the DMPK 
Numerical simulations EBI. 



equation. Numerical simulations ||46|J47| | showed that the correlations between a pair of a m 
are rather weak. So it might be that the correlations vanish in the thermodynamic limit 
leading to a diagonal form of A mn - A perturbative calculation of the limiting distribution 
for strong forward scattering will be published in a subsequent paper |35| . It will shed some 
light on these questions. 
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APPENDIX A: DERIVATION OF THE LANGEVIN EQUATIONS 

The Langevin equations for T and u which describe the stochastic evolution of the matrix 

TT^ = M-( Mn Ml2 ) 
11 ~ m ~ \ M 21 M 22 J 

( it(cosh T)u^ n(sinh r)n T \ , . . 

= ( u*(sinhr)ut u*(coshT)u T I ^ > 



were already given by Dorokhov in his pioneering work |32|| . Since he did not derive them 
explicitly we derive them in this appendix. The multiplicative nature (|I5|) of the transfer 
matrix implies that 

M(L + 5L) = T(L + 5L, L)M(L)T\L + 8L, L) (A2) 

if a short segment of length 5L is added to a sample of length L. The change AM = 
M(L + 5L) — M(L) of M induces the changes AT and Au which can be calculated by 
perturbation theory. The eigenvalues of the hermitean matrix M 11 = ncoshT are 
coshT m . The corresponding eigenvector to coshr m is the m th column vector M m of u. The 
change An may be expanded into these eigenvectors 

JV 

Au m = c mk u k . (A3) 

k=l 

Non-degenerate first-order perturbation theory then yields 
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Ar m = {tttA y" + 0((AM 11 ) 2 ) (A4) 



and 

{utAM%} ™ +0((AM-) 2 ) (A5) 



cosh 1 m — cosh 1 rt 
for n m. The expansion coefficient c mm can be calculated from 

AM 12 = Au(F + AF)n T + u(F + AF)Au T + uAFu T + Au(F + AF)Au T (A6) 

where .F = sinhT. Equation (|A3|) implies that m^Au = c T . Multiplying Eq. ( TA~6| ) with 
from the left and with u* from the right thus yields 

v)AM l2 u* = c T (F + AF) + (F + AF)c + AF + c T (F + AF)c . (A7) 

Hence 

^ = {.WU ^i^H. + Q( ( AM")») . (A8) 

smh 1 m 

Inserting the expansion (|21| ) of T(L + <5L, L) into powers of 5L into Eq. ( |A2|) gives 

AM ij = (y k M kj + M ik ~/ kj f ) + 0(5L 2 ). (A9) 
Collecting results and taking the limit <5L — > finally leads to the Langevin equations 

, ^ f E kn sinh T n + E* kn sinh T k 



k k^n \ 



dL h m t^v V cosn — cosn Tfc 



Umk 



+ Q coth r„(£ n „ - £*J^ (A10) 

where E = u + ~f 12 u*. Note that the symmetries y 11 ^ = — 7 11 and ^ 12T = -y 12 imply that 
d{wu)) / dL\ uu \ =1 = which ensures that unitarity is conserved. 



APPENDIX B: THE COEFFICIENTS OF THE FP OPERATOR 

The coefficients of the FP operator can be calculated by iterative integration of the 
Langevin equations (^) and subsequent averaging over the disorder |57| . Integration of the 
Langevin equation yields 

T m {x) = T m {L) + J* dx\E mm {x') + E* mm (x')) 

(L) + jf dx'{j u (x')u(x') + u{x')P{x')} mn . (Bl) 
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The matrix elements of u(x') and u*(x') which appear in the integrands can again be 
expressed by the second of the Eqs. (|B1|) . Iterating this procedure leads to an increasing 
number of terms with polynomials of 7^ of increasing degree r. For the Gaussian white noise 
model (|33|) only polynomials of even degree have non-zero disorder averages. The integration 
over the 5-functions of the disorder averaged polynomials then leads to terms of order 5L r ^ 2 . 
Hence, only polynomials of degree two contribute to the limit [- • •] = lim^^of ■ -)sl/SL. 
Collecting these polynomials leads to the following result for the coefficients of the FP 
operator 



[Ar, 
[Ar m Ar T 

[AT m Aii m / n 
[Ait mr 



[Au mn Au 



[Aw mn Ait m , n /] 



where 



/•L+SL r-L+dL 

[ j L dx J l dx>V 2 (s)7 U *(x')u* + £(x)V*(x') + c.c.} r 

rL+8L pL+SL 

[ J l dx J l dx'{£(x) + £*(x)} mm {£(x') + £*(x')} nn ] 

(■L+5L i-L+SL 

[ J dx J dx'{£(x) + £*(x)} mm { 7 n (s')n + uV{x'))} m , 



-L+8L 



- [ / dx dx'{-f u (x) 7 n (x')u + 2 7 n (x)uV(x') 

2 JL JL 

+uV(x')V(x)} mn 

+ E U ™A™ ~ 3){U^ 1X V)7 12 (SX 



rL+8L 



+u T 7 12 (x) 7 u *(x'K + V\x')£{x) + £{x)V*{x')} jn 
-J2 u ™j°U - n){u^j uj (x')j l2 (x)u* 



u^ 12 (x)j 11 *(x')u* + V\x')£{x) + £{x)V*{x')}* n 



11 * / /\ * 



/L+8L rlj+dlj 
dx J dx'{j u (x)u + uV(x))} r 

{7 (x')u + uV[x 

) ) jra'n'] 



rL+8L 



r-L+8L rL+OL 

I dx dx'{j n (x)u + uV(x))} r 

JL JL 

{ 7 11 *(x / )n* + n*P*(x / ))}w] 



r-L+SL 



(B2) 



£{x) = u^ 12 (x)u* 
V m n{x) = 9{n - m)S mn (x) - 9(m - n)£^ n (x) 



and integrals of the type f L L+5L dx J^dx'S(x — x')f(x') have been replaced by J 
^ +8L dx'5{x-x')f{x')/2. 



(B3) 
t +SL dx 
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APPENDIX C: THE INVARIANT MEASURE OF THE UNITARY GROUP 



The invariant measure dfi(u) of the unitary group is invariant under multiplication with 
an arbitrary element Uq of the group from the left and the right 

dfi(u) = dfi(u u) = dfi(uu ). (CI) 

As claimed in section |TTT] we show in this appendix that the invariant measure has the form 

dfi(u) = V-\N)5(l,uu^ J] dx mn dy mn , (C2) 

m,n 

where u mn = x mn + iy mn , V(iV) is the volume of the unitary group, and the ^-function 
6(1, uu^) has been defined in Eq. (£T). Since the unitary group is compact, invariance 



under multiplication from the right implies invariance under multiplication from the left (cf. 
p. 316 in |[59[ ). Therefore it is sufficient to show the right invariance. Let us write u' = uuq. 
Then 



dp(u') = V-\N)6(1, u'u t) J] dx' mn dy' m 



m,n 



V~ 1 (N)6(l,uu i ) 



det 



d(x mn , Umn) 



dx mn dy mn . (C3) 



m,n 



Thus dfi(u) is right invariant, if the absolute value of the Jacobi determinant of the linear 
transformation u' = uuq equals one. This transformation is equivalent to 

/ x \ = / 1 N ® xl -ljv ® Vq \ ( x\ 
\ y ) I ljv ® yl l N ®xl ) \y ) 



d(x, y) 



x 



(C4) 



9(iB,y) V W / 

where u = x + iy has been written in the vector form 

x T y T ) = ( a^ii ^12 • • • a;iiv ^21 • • • x NN y n • • • 2/jvjv ) • (C5) 



The multiplication rule (A ® B)(C ® D) = (AC) ® (BD) for cross products of matrices 
and the unitarity of it imply that d(x , y )/d(x, y) is an orthogonal matrix. Therefore the 
absolute value of the Jacobi determinant is one and the measure is right invariant. 

APPENDIX D: ALTERNATIVE DERIVATION OF THE FP EQUATION 



The FP equation ( p9|) can be derived in an alternative way yielding immediately the form 
in which the derivatives act directly on the probability distribution and not on the coefficients 
of the FP operator. This implies useful relations between these two forms which are difficult 
to obtain by direct differentiation of the coefficients if the result is not known in advance. 
Therefore we describe the alternative derivation in this appendix. It is a generalization of 
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the technique which was employed by Mello and coworkers to derive the FP equation of 
their isotropic model [jS0j , [l0|l . 

Assume that the probability distribution of M for a system of length L is p(L;M). 
Then add a statistically independent segment of length SL to the system. The distribution 
of the transfer matrix Tsl of the segment is denoted by w(L, 6L; Tsl)- Averaging a function 
f(M) over the disorder of the whole system yields 

(f(M)) L+SL = J 'dp (M) p{L + 5L- M) f(M) 

= J J dp (M 1 ) du (T SL ) P(L; M ') w(L, 5L; T 5L )f(T 5L M' T\ L ) (Dl) 

where M = T SL M' T\ L and duo (T), dp (M) are measures on the matrix spaces of T and 
M. If the measure dp (M) is chosen to be invariant under the transformation T M T\ for 
any transfer matrix T one finds 

(f(M)) L+SL = JJdp (M) duo (T 6L ) p(L; T^M{T\ L Y l )w{L, 5L; T sl ) f(M). (D2) 

A similar line of reasoning as in JlO| yields that the invariant measure has the form 

N 

dp (M) = J(T) J] dT m dp(u) (D3) 

m=l 

where J(T) = Y[ m<n \ coshr, m — coshT n | rimS^hrm and dp{u) is the invariant measure on 
the unitary group. Comparing Eq. ( pi|) with Eq. ( P2| ) immediately shows that 

p(L + 6L; M) = J duo (Tsl) P(L; T^M (T^)' 1 ) w(L, 5L; T sl ) (D4) 

or equivalently 

p(L + 5L; M) = (p(L; Tf L M (T\ l Y 1 )) sl (D5) 

since the average of p(L; TjlM (Tj L ) _1 ) with the probability w(L, 5L; TsL)duo(TsL) is just 
the average (• ■ -)sl over the disorder of the segment. For a fixed value of T the distribution 
p(L + SL; M) may be expanded into the irreducible representations of the unitary group. 
Since the irreducible representations are polynomials in u mn and u* mn |6T| , |62| , p(L + 5L; M) 
can be analytically continued to arbitrary complex matrix elements u mn . This justifies to 
write Eq. ( p5| ) in the form 

p(L + 5L- T, u, u*) = (p(L; T + AT, u + Aw, u* + Au*)) SL , (D6) 



where AT, An and An* are the changes of the parameters of M which are induced by the 
mTjlM 
in power 

dp ( 1 



transformation Tj^M (TJ l ) _1 . A Taylor expansion of the left side in powers of 5L and of 
the right side in powers of Ar m , Au mn and Au* yields 



m l2 [Ar m Ar n ]9r m «9r„ + [Ar m ]dr n + [Ar m Au m , n ,}dr m d Um>nl 
+[Ar m Au* m , n ,)dr m d u * m/ni + [Au mn )d Umn + [Au* mn ]d u * mn 

+ [Au mn Au* m , n ,]d Umn d u * nini + -[Au mn Au mln ]d Umn d Um , n , 

+ll^<n^<>n>}du* mn du llnl }p (D7) 
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in the limit 5L — > 0. The coefficients [•••]= limsL->o(- • -)sl/8L could be determined by the 
calculation of Ar m , Au mn and At^ n with perturbation theory and subsequent averaging 
over the disorder. This would lead to expressions which involve [e^J and [£^ n £^y n /] where 

T(L + 8L,L)- x = l+[ e " £ Z\, (D8) 



£ £ 



-11. 



e XJ -=£ n ', and £ 12 = —£ 12T . In order to include all the terms which contribute to the 
coefficients one had to go to the second order of the perturbation theory which is quite 
involved. We proceed in a different way, instead. For the Gaussian white noise model 
one finds 



[4J = [e 



»J 1 

mnl 



L mn m'n 



[ £ mn £ rn'n>]- (D9) 



Hence, the coefficients may be as well evaluated with the changes Ar m , Au mn and Au* mn 
which are induced by the transformation T $lMT\ l . These changes can be obtained by 
iterative integration of the Langevin equations ( |A10| ) similar to that described in appendix 
|B| for the simpler case of large system lengths. 

Since dp(M) / J(T)=Y[ m dT m dfi(u) is the same measure which was used for the probabil- 
ity distribution of the FP equation (|39|), one expects that the distribution J(T)p(L; T, u, u*) 
obeys this FP equation for large system lengths. In fact, the transformation of Eq. ( [Q7| j to 
this distribution and the neglect of the exponential small contributions to terms of the form 
(Pff) leads to Eq. (^) where the operators A, B and C appear in the form 

A mn = - [Ar m Ar n ] 

B m = [AT m ] - n [Ar n Ar m ] 

+ [Ar m Au w ] d Um , n , + [Ar m A< n , n ,] d u * mfn , 

C = ^mn [Ar m Ar n ] - m [AT m ] - m [Ar m A?i m / n /] d Um , n , 
-m [Ar m A^,„,] d u * mini + [Au mn } d Umn + [Au* mn ] d u * mn 

+ [Au mn Au* m , n ,\ d Umn d u * mini + - [Au mri Au mlril ] d Umn d Um , n , 
+hAu* mn Au* m ,Jd u * mn d u * (D10) 

The equivalence with the form of A, B and C in Eq. (f40D can be shown by a lengthy but 
straightforward calculation which exploits only the symmetries 7 11 ' = — 7 11 and ^ 12T = 
7 12 - 
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